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Finding cost-optimal solutions towards nearly-zero-energy buildings (nZEBs) in accordance with Euro¬ 
pean energy performance of buildings directive (EPBD-recast 2010) is a challenging task. It requires 
exploring a huge number of possible combinations of energy-saving measures (ESMs) and energy-supply 
systems including renewable energy sources (RESs), under a comparative framework methodology. The 
current study introduces efficient, transparent, and time-saving simulation-based optimization method 
for such explorations. The method is applied to find the cost-optimal and nZEB energy performance levels 
for a study case of a single-family house in Finland. Different options of building-envelope parameters, 
heat-recovery units, and heating/cooling systems as well as various sizes of thermal and photovoltaic 
solar systems are explored as design options via three-stage optimization. The resulted economic and 
environmental trade-offs show that primary energy consumption >93 and <103kWh/m 2 a is a cost- 
optimal energy performance level. It is economically feasible to achieve nZEB with 70 kWh/m 2 a. However, 
incentives (e.g., energy credits) are required to reach lower-environmental-impact houses. Investing in 
low-operating-cost environmentally friendly heating system (e.g. ground source heat pump) is a key 
element for optimal solutions. The optimal implementation of ESMs and RES depends significantly on 
the installed heating/cooling system and the escalation rate of the energy price. 

© 2012 Elsevier B.V. All rights reserved. 


1. Introduction 

Buildings are responsible for 40% of energy consumption and 
36% of the EU’s C0 2 emissions. Energy performance of buildings is 
a key element to achieve the EU climate and energy objectives, 
namely a 20% reduction of the greenhouse gases emissions and 
20% of primary energy savings by 2020. Improving the energy per¬ 
formance of buildings is a cost-effective way of fighting against 
climate change and improving energy security [1-3]. From the 
economic point-of-view, the best environmental solutions cannot 
be guaranteed by regulations which depend mainly on building 
envelope requirements [4-7]. The analysis of energy efficiency and 
cost optimality should consider the energy sources and building 
systems. Seeking for cost-optimal high-energy-performance build¬ 
ings, the EPBD-recast 2010 [8] requests the Member States to avoid 
establishing rules whereby a measure on the building envelope is 
always applied first and only then a measure on a building system 
is allowed. According to the EPBD recast, the minimum energy per¬ 
formance requirements should be set with a view of achieving cost 
optimal levels for buildings, building units and building elements. 


* Corresponding author. Tel.: +358458747044; fax: +358 9 470 23418. 
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Fligher energy performance buildings, like nZEBs, should also be 
economically feasible. Finding cost-optimal minimum energy per¬ 
formance requirements and nZEB solutions is an arduous task. The 
task requires exploring a huge number of design solutions (com¬ 
binations of energy saving measures and energy supply systems) 
under a comparative framework methodology. The EPBD recast 
required the EU Commission to establish a comprehensive method¬ 
ology by 30th June 2011.The methodology is demonstrated by the 
Buildings Performance Institute Europe (BPIE) [9] and the Euro¬ 
pean Council for an Energy-Efficient Economy (ECEEE) [10]. The 
current study introduces a multi-stage optimization method for 
exploring wide spaces of building and system integrated solutions, 
transparently and efficiently. The method is designed to reduce the 
exploration and analysis efforts needed to find optimal solutions for 
the new EU-buildings in line with the EPBD framework method¬ 
ology. The method uses simulation-based optimization approach 
hence it should be suitable for most of new buildings in Europe 
where the heating is the major demand for thermal comfort. 

2.1. The EPBD recast comparative framework methodology 

The EPBD recast comparative framework methodology was 
established for calculating cost-optimal levels of minimum energy 
performance requirements for buildings and building elements 
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Nomenclature 

AHU 

air handling unit 

COPcu 

coefficient of performance of the cooling unit 

dEle 

difference between electricity demand and produc¬ 
tion [l<Wh/m 2 a] 

dLCC 

difference in life-cycle costs between any solution 
and the reference design [€/m 2 ] 

dPV 

difference in present value between any solution 
and the reference design [kWh/m 2 a] 

dQDHW 

domestic hot water energy saving by the solar ther¬ 
mal collector [kWh/m 2 a] 

dh 27 

degree hours over 27 °C indoor temperature [°C h] 

DH 

district heating 

DHW 

domestic hot water 

DHW ele 

the electrical portion of domestic hot water 
[kWh/m 2 a] 

Ehv 

the electricity consumption of the HVAC systems 
(fans and pumps) [l<Wh/m 2 a] 

Ela 

the electricity consumption of lighting and appli¬ 
ances [kWh/m 2 a] 

Ele 

electricity consumption [kWh/m 2 a] 

EH 

electrical heating 

EPBD 

energy performance of buildings directive 

ESMs 

energy saving measures 

GSHP 

ground source heat pump 

IC 

investment cost [€/m 2 ] 

LCA 

life-cycle assessment 

NSGAJI 

elitist non-dominated sorting genetic algorithm 

nZEBs 

nearly zero energy buildings 

OB 

oil boiler 

OC 

operating cost of energy [€/m 2 ] 

PEC 

primary energy consumption [kWh/m 2 a] 

PV 

photovoltaic 

PV e 

the useful electricity produced by photovoltaic sys¬ 
tem [kWh/m 2 a] 

PW 

present worth 

Qc 

space cooling energy demands [l<Wh/m 2 a] 

Qh 

space-heating energy demands [kWh/m 2 a] 

RC 

replacement cost [€/m 2 ] 

RES 

renewable energy sources 

SH 

space heating [l<Wh/m 2 a] 

SH e i e 

the electrical portion of space heating [l<Wh/m 2 a] 

SHGC 

solar heat gain coefficient 

SWSC 

short wave shading coefficient 

S-factor 

solar gain factor 

T-value 

solar transmittance 

H-value 

heat-transfer coefficient [W/m 2 I<] 

?7shs 

efficiency of the space-heating system [%] 

?7dhws 

efficiency of the domestic hot water system [%] 

delist 

distribution efficiency of the heating system [%] 


[9,10]. The methodology requires comparing the global costs (addi¬ 
tional investments, replacement costs, energy costs, etc.) and the 
delivered primary energy of combinations of compatible energy 
efficiency and energy supply measures (packages of measures). 
The packages should range from those in compliance with the cur¬ 
rent regulations to combinations that realize nZEBs. The packages 
should also include various options for renewable energy sources 
(RES) generation. Fig. 1 shows a predicted cost-optimal curve that 
will be found when assessing all combinations of commonly used 
and advanced measures. The lowest part of the curve represents the 
economic optimum for the combinations of measures. The min¬ 
imum energy performance requirements are represented by the 


Global 

costs 

[€/m 2 ] 


, 4 



Distance to target 
(new buildings, 2021 


Primary energy 

consumption 

[kWh/m 2 ] 


Fig. 1 . Cost optimal curve and distance to 2021 target [9]. 


economic optimum point that delivers the lowest cost for the end- 
user and/or for the company or society. The part of the curve to the 
right of the economic optimum represents solutions that underper¬ 
form in both aspects (environmental and financial). The left part of 
the curve, starting from the economic optimum point, represents 
the cost-optimal energy-performance levels for low and nearly- 
zero-energy buildings. The figure shows also the distance to the 
EU-2021 target (nZEB) for new buildings. 


2.2. The aim of the current study 

One of the main challenges of the EBPD-2010 calculation 
methodology is to ensure that, on the one hand, all measures with 
a possible impact on the primary or final energy use of a building 
are considered, whilst, on the other hand, the calculation exercise 
remains manageable and proportionate [11]. Applying few options 
for several variants could offer millions of design solutions. In order 
to limit the number of solutions, the guideline of the EBPD draft 
[11] proposes to address a matrix of energy efficiency packages, 
which rules out mutually exclusive technologies. For instance, a 
heat pump for space heating (SH) does not have to be assessed in 
combination with a high efficiency boiler for space heating as the 
options are mutually exclusive and do not complement each other. 
The possible energy efficiency measures and measures based on 
RES (and packages/variants thereof) can be presented in a matrix 
and unfeasible combinations eliminated. 

The elimination approach cannot guarantee global cost-optimal 
solutions because it explores only some of the available combi¬ 
nations of design options. Furthermore, considerable effort and 
experience are needed to make correct eliminations. To estab¬ 
lish a comprehensive overview, all compatible combinations of 
commonly used and advanced measures should be considered. 
Stochastic methods are promising, and can be used to investigate a 
huge number of combinations. However, they should be employed 
under a suitable scheme. The aim of this study is to introduce 
a suitable optimization scheme/method which provides efficient, 
transparent, and time-saving exploration. 

• Efficient exploration is performed by using a combination of a 
modified elitist multi-objective non-nominated sorting genetic 
algorithm and detailed simulation programs. 
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• Transparent exploration is presented via three optimization 
stages showing the effect of the design-variable combinations on 
the objective and constraint functions. 

• Time-saving exploration is achieved via speeding up the explo¬ 
ration by avoiding the unfeasible design-variable combinations 
and using pre-simulated results instead of running time- 
consuming simulation (when possible). 

1.3. Relevant studies 

The EBPD recast-2010 is still a quite new directive. Therefore, 
it has not been widely complied with. The economic environment 
scope of the EBPD recast is based on the European Standard EN 
15459 [12]. BS ISO 15686-5 [13] is the first international standard 
for property life-cycle costing. The EN 15459 and ISO 15686-5 were 
complied with in [14] and [15], respectively. Many of building 
studies, e.g., [14-18] investigate the economic and environment 
point-of-views. However, the majority address a limited number 
of design solutions. For instance, Georges et al. [14] investigated 
a single-family dwelling in Belgium, addressing combinations of 
sixteen heating systems and five building designs. Marszal and 
Heiselberg [15] found the minimum life-cycle cost for a multi¬ 
storey residential net-zero-energy building in Denmark, addressing 
three levels of energy demand and three alternatives of energy 
supply systems. In Finland, the first house designed for mini¬ 
mal energy consumption was built in the early 1990s, and its 
energy consumption was monitored for three years. On the basis 
of this monitoring, VTT (Technical Research Centre of Finland) 
estimated that additional investment in energy efficiency pays 
itself back in 5-6 years and yields a return of 15-20% [19]. VTT’s 
energy renovation technologies project [20] studied the profit¬ 
ability of energy renovation measures (structural improvements, 
heat supply systems, ventilation systems, lighting, electrical appli¬ 
ances, solar shading and cooling) for three types of buildings. The 
project applied the simple payback time method for its calculation. 
Pylsy and Kalema [21] made a life-cycle cost sensitivity analysis, 
addressing four building insulation levels, four building tightness 
levels, three ventilation-heat recovery types, and nine heating sys¬ 
tems to find cost-effective concepts for low-energy houses. The 
analysis found that the improvement of the thermal insulation 
of the building envelope is the most effective way to reduce the 
space-heating energy need. However, the lowest building insula¬ 
tion level is selected as a cost-optimal solution when ground-source 
heat pump (GSHP) is selected for heating. Hasan et al. [22] com¬ 
bined simulation and optimization to minimize the life-cycle cost 
(LCC) of a single-family detached house. The study investigated a 
wide range of wall, roof, and floor insulation levels, two types of 
windows, and two types of ventilation-heat recovery units. How¬ 
ever, no heating alternatives were addressed. Alanne et al. [23] 
considered the selection of a residential energy-supply system as 
a multi-criteria decision-making problem involving both financial 
and environmental issues. The study analyzed the competitiveness 
of micro-CHP as 1 of 10 alternative heating systems for a Finnish 
single-family house. The analysis showed that the micro-CHP is 
a reasonable alternative to traditional systems, particularly from 
the environmental point of view. Alanne et al. [24] searched for 
optimized strategies for the integration of a Stirling engine-based 
micro-cogeneration system in residential buildings by comparing 
the performance of various system configurations and operational 
strategies with that of a reference system, i.e., hydronic heating and 
a low temperature gas boiler in standard and passive house con¬ 
structions located in different climates. Saari et al. [25] studied eight 
different energy-saving design concepts and three heating modes 
for a typical new Finnish detached house. The study showed that 
the payback time of the heating system depends on the real interest 
rate and the building construction. 
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Fig. 2. Plan of the two-storey house studied (143 sq m). 

The above Finnish studies focused mainly on the heating energy. 
This is not consistent with the EBPD recast 2010, which focuses on 
the total energy consumption (including the energy consumption of 
heating, cooling, lighting, appliances, and auxiliary equipment) and 
its financial and environmental impacts. Generally, the reviewed 
studies either are inconsistent with the EPBD-recast 2010 or cannot 
guarantee global cost-optimal solutions because they investigate 
a limited number of design solutions. As mentioned above, the 
current study introduces efficient, transparent, and time-saving 
exploration for wide solution spaces. The introduced method is 
applied to find the cost-optimal and nZEB energy performance lev¬ 
els for a single-family house in Finland in line with the EPBD recast. 

2. Methodology 

2.1. Defining the building and its conceptual design 

Energy efficiency of houses can strongly influence overall energy 
use in the building sector. In Finland, residential buildings com¬ 
prise 85% of the building stock and occupy 64% of the building area 
[25]. A single-family house is chosen as a residential building case 
study (Fig. 2) because about 70% of the Finnish population live in 
houses [26]. The house, in Helsinki (60° N, 25° E), consists of two 
storeys. The floor area of the house is 143 sq m. The internal height 
of each storey is 2.5 m. The glassing area represents about 15% of the 
floor area. The two storeys are connected by a staircase and have 
operable windows for natural cooling in summer. To avoid sum¬ 
mer overheating risk, different types of shading and a mechanical 
cooling unit are offered as design options. Building tightness and 
insulation levels could also control the risk. All the building enve¬ 
lope and system options are described in Section 2.3. According to 
the Finnish Building code D5-2012 [27], the construction type of 
the house is Medium-heavy I. The house has interior effective ther¬ 
mal capacity of about 21 kWh/K. The house is ventilated by one 
air-handling unit (AHU), which supplies fresh air to the bedrooms 
and living room and draws the exhaust air from the bathrooms and 


Living Room 


Stair 


WC 

Kitchen 

Bed R. 




\ 


Stair 


Small window for 
summer cooling 
<0.9 m 2 ) 






































192 


M. Hamdy et al. / Energy and Buildings 56 (2013) 189-203 


Table 1 

Insulation thickness (design variables no. 1, 2 and 3). 


Insulation of... 

Material 

Range of insulation 
thickness (m) 

No. options using a 
uniform step 
(0.02 m) 

U-value of the 
construction 

(W/m 2 K) 

Investment cost 

External wall 

Mineral wool 

From 0.185 to 0.48 

16 

From 0.17 to 0.07 

62.72 

Roof 

Blow-in wool 

From 0.41 to 0.55 

8 

From 0.09 to 0.07 

36.2 

Floor 

Polyure-thane 

From 0.2 to 0.44 

13 

From 0.17 to 0.080 

111.4 


Costs are taken from [7], [36] and updated to 2011 ones by using 3.8% inflation rate [37] considering the rise in the material costs. 


Table 2 

Building tightness levels (design variable no. 4). 


Levels 

Specification n 50 (1/h) 

Price for labour cost (€/m 2 ) 

1 

2 

0 

2 

1 

11 

3 

0.5 

22 


Costs are taken from [7] and updated to 2011 ones by using 2.1% inflation rate [37] 
considering the rise in the labour costs. 


the kitchen. The AHU heating coil keeps the supply air tempera¬ 
ture at 18 °C when the incoming outdoor air temperature is lower 
than this temperature. The average exhaust-air flow from the whole 
house is equal to 0.65 air change per hour [1/h], which is higher 
than the minimum requirements (0.5 1/h) of the Finnish Building 
code D2-2012 [28]. The internal gains due to people, lighting, and 
electric appliances are 10.3, 7.8, and 17.81<Wh/m 2 a, respectively. 
The values are in accordance with the D5-2012 [27]. The dynamic 
changes of domestic hot water, lighting, and appliance energy are 
considered by using profiles based on the typical Finnish life style 
[29]. The domestic hot water tank has the size of 3001. The supply- 
water temperature to the space heating is controlled as a function of 
the outdoor air temperature. The electrical-energy consumptions of 
the auxiliary equipment (fans and pumps) are calculated, assuming 
0.625 l<W/(m 3 /s) specific fan power and 151<Pa nominal pumping 
pressure head. The fans and pumps have efficiencies of 60% and 
50%, respectively. Both consume about 7.5 kWh/m 2 a of electricity. 
The dynamic thermal performance of the house and systems are 
simulated by IDA ICE [30] and IDA ESBO [31]. The simulation uses 
reference year weather data (Vantaa_TRY2012) which is developed 
recently for energy calculation purposes in Finland [32]. 



Ground source heat 
Exhaust air heat pump 

■ Electric heating 

■ Oil fired heating 
District heating 

■ Pellet heating 

■ Other forms of wood 
heating 


Fig. 3. The market share of heating systems in new single-family houses in Finland 
[34]. 


are heated by electricity (Fig. 3). Because of the mild climate in sum¬ 
mer, most of the Finnish houses do not have mechanical cooling. 
In the current study, the reference design is served by electri¬ 
cal radiators for space heating; electrical boiler for domestic hot 
water; and mechanical ventilation unit with 60% heat recovery effi¬ 
ciency. Neither mechanical cooling nor on-site renewable systems 
are implemented. The calculation by the IDA ICE [30] building simu¬ 
lation program indicates that the required space-heating energy for 
the reference case is 461<Wh/m 2 a. The annual domestic hot water 
demand is 32 l<Wh/m 2 a. The total electricity demand of the house 
is 126.5 kWh/m 2 a, including fans, pumps, appliances, lighting, and 
heating demands. 


2.2. Defining the building reference design 

In line with the EPBD comparative framework, the economical 
and environmental viability of the design solutions is assessed rel¬ 
ative to a reference design. The reference design is with D-values 
0.17,0.09,0.17, and 1 W/m 2 K for external walls, roof, ground floor, 
and window, respectively. The values are based on the National 
Building Code of Finland, C3-2010 [33]. The house has a building 
tightness (n 50 ) of 2 1/h, where n 50 is the number of air changes 
per hour equivalent to an air-leakage rate with a 50 Pa pressure 
difference between indoor and outdoor. The window type has a 
solar heat gain coefficient (SHGC) and a solar transmittance (T- 
value) of 46 and 34%, respectively. To achieve acceptable level of 
summer overheating, the windows are shaded by blinds, with hor¬ 
izontal laths, between the inner panes. According to the market 
share of heating systems [34], 42% of the small houses in Finland 


2.3. Defining the design variables 

The design variables are selected to cover packages of meas¬ 
ures ranging from compliance with the requirements of the current 
national building code C3-2010 [33] to combinations that realize 
nZEBs (e.g., passive house D-values [35], photovoltaic, and solar 
thermal collectors). The variables include a number of external 
wall, roof, and floor insulation thicknesses (Table 1), three build¬ 
ing tightness levels (Table 2), three window types (Table 3), four 
shading options (Table 4), three heat recovery units (Table 5), two 
cooling options (Table 6), four heating systems (Table 7), and dif¬ 
ferent sizes of on-site solar systems (Table 8). The costs of the 
ESMs (insulation, building tightness, window type, and heat recov¬ 
ery unit) are from previous studies [7,36]. The costs are updated 
here by using building-cost indices and inflation rate from [37]. 
The technical properties of window and shading options are taken 


Table 3 

Window types (design variable no. 5). 


Type 

U-value (W/m 2 K) 

T-value 

SHGC 

Cost (€/ m 2 ) 

Description 

1 

1 

0.34 

0.46 

228.4 

Triple Laminated glass Wood aluminum frame (Argon gas) 

2 

0.85 

0.29 

0.42 

267.4 

Triple Laminated glass Wood aluminum Frame (Argon gas) 

3 

1.1 

0.28 

0.38 

234 

Quadruple Laminated Wood aluminum frame(Argon gas) 


Costs are taken from [7] and updated to 2011 ones by using 3.8% inflation rate [37] considering the rise in the material costs. Technical description [38]. 











Table 4 

Shading options (design variable no. 6). 
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Type 

S-factor 

SWSC 

Description 

Cost (€/ m 2 ) 

1 

0.14 

0.09 

External blind, horizontal laths 

200 

2 

0.33 

0.12 

Blind between the outer panes, horizontal laths 

60 

3 

0.53 

0.14 

Blind between the inner panes, horizontal laths 

50 

4 

0.65 

0.16 

Internal blind, horizontal laths 

25 


Costs are based on the Finnish market 2011 prices. Technical description: IDA ICE [30]. 

The windows shading decreases SHGC by (1-S factor)*! 00% and T-value by (1 - SWSC) x 100%. 


Table 5 

Heat recovery unit (design variable no. 7). 


Type 

Annual efficiency (%) 

Description 

Cost (€/ m 2 ) 

1 

60 

Cross-flow heat exchanger 

3533 

2 

70 

Counter-flow heat exchanger 

3835 

3 

80 

Regenerative heat exchanger 

4138 


Costs are taken from [7] and updated to 2011 ones by using 2% inflation rate consid¬ 
ering the rise in the system costs. 


from [30,38], respectively. The installation costs of the heating and 
cooling systems are assumed based on the Finnish market prices 
[39-44]. The mechanical cooling option (Table 6) with a set-point 
of 25 °C is implemented for building designs which cannot provide 
comfort conditions [45] in summer by shading (Table 4) and/or nat¬ 
ural ventilation via window opening. The energy prices are taken 
from [46,47]. 

Since solar energy is one of the most favorable RESs in Europe 
[48], the current study offers solar thermal and photovoltaic sys¬ 
tems to reduce the DHW and electricity demands. An economic 
feasibility study of a net-zero-energy house in the cold climate of 
Montreal [49] concluded that flat-plate solar collectors are the best 
financial option regardless of how many collectors are installed. 


Flat panels are the most durable type of collectors, and they also 
have the best performance for systems designed for temperatures 
within 100 degrees. Based on this conclusion, the current study 
investigates only the flat plat collector as a cost-effective solar ther¬ 
mal option. The installation cost of a solar thermal kit includes the 
costs of a flat plate collector, a storage tank with heat exchang¬ 
ers, a circulating pump, piping and insulation. The 1000 € cost of a 
regular storage tank (3001), to be installed on a non-solar thermal 
system, is deducted from the total cost. The price of a photovoltaic 
module and inverter are based on the average prices in Europe (2.37 
and 0.53 €/W, respectively) [50]. By adding the installation cost, the 
total cost reaches 6€/W. 

2.4. Defining the objective functions 

In order to find the cost optimal curve (Fig. 1 ), a multi-objective 
optimization problem is defined with two objective functions: 

Min{/!(x),/ 2 (S)} x = [x u x 2 , . x m ] 

where 

f \: primary energy consumption PEC. 


Table 6 

Cooling options (design variable no. 8). 


Options 

Specifications 

COPcu 

Cost (€/unit) 

1 

No cooling 

- 

0 

2 

Small unit 2.5 kW cooling capacity 
with rotary piston compressor, 
refrigerant R407C Eco 

3.2 

850 


Costs are based on the Finnish market 2011 prices [44]. 


Table 7 

Primary heating systems (design variable no. 9). 


System 

Investment cost (€) 

Service cost 
(€/a) 

Subscrip-tion 
fee (€/a) 

Energy price 
(c/kWh) 

OSHS (%) 

? 7 dhws (%) 

^7dist (^) 

Energy factor (F) 

Direct electricity with 
electrical radiators (EH) 

50 kW p + 2700 

30 

83 

13.5 10.9 a 

100 

88 

94 

1.7 

Oil boiler with water 
radiators (OB) 

286 kW p +7143 

135 

83 

6.12 

81 

81 

87 

1 

District heating with water 
radiators (DH) 

50.5 kW p +9050 

40 

404 b 

6.5 

94 

94 

87 

0.7 

GSHP with floor heating 
(GSHP) 

592.5 kW p +12155 

145 

83 

13.5 10.9 a 

300 

250 

84 

1.7 


Energy prices are based on the Statistics Finland [46] and Helsingin Energia [47]. 

Costs are taken from [39-42] and updated to 2011 ones by using 2% inflation rate. 
a The price of day electricity (13.5 c/kW) on weekdays, from Monday to Friday, 7 am-8 pm. The price of night-time electricity (10.9 c/kW) at other times [41 ]. 
b Beside the 83 € annual fee of electrical connection, 321 €is added for district heating connection. 


Table 8 

On-site RES supplementary systems (design variables no. 10 and 11). 


Options 

Solar system 

Total cost (€) 

Size (sqm) 

1 

Flat plate solar thermal system 

380 Acollector + 3200 

From 0 to 30 m 2 

2 

Multicrystalline silicon Photovoltaic (PV) system 

6€/W 

From 0 to 70 m 2 


Costs are based on the Finnish market 2011 prices [43,50]. 
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Fig. 4. The annual domestic hot-water energy saving (c/Qdhw) using different areas 
of flat-plate solar thermal collector. 
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f 2 : difference in life-cycle cost between any design option and 
the reference design dLCC. 

x: combination of the design-variables (xi,x 2 ,.. .x m ). 
m: number of the design variables. 

2.4.2. Primary energy consumption (PEC, first objective function) 

In order to achieve environmental solutions, the PEC is mini¬ 
mized as a first objective function. Primary energy is the favorable 
metric, according to the latest EPBD recast [8 ]. The PEC considers the 
energy use of the building (heating, cooling, ventilation, lighting, 
pumps and fans, other technical service systems, DHW, cooking, 
appliances, lighting) and the energy-saving by RES generations. The 
PEC is calculated by using non-renewable primary energy factors F 
according to the energy source (Table 7) 

PEC = F SHdeijyej-gd + F DHW^ e ij verec j + F Ele^ e ii verec j (1) 

cu _ (Qh/^dist) 

^^delivered — ~ 

4shs 


DHWdeijygj-ed — 


(Qdhw - dQpHw) 

4DHWS 


dEle(f) = + SHeie(t) + DHWeie(t) + Ehv(t) + E la (t) - 

Cu vinverter 


t=8761 

Ele de iivered = ^ max {dEle(f), o} 

t=l 

Eq. (1) divides the energy demands (Qh, Qdhw. Qc) by the annual 
efficiencies to calculate the delivered ones. According to the heat¬ 
ing system applications (space heating and the domestic hot water 
applications), two efficiencies (rjsHS, and 4 dhws) are considered as 
being consistent with the Finnish regulation [27]. Based on the 
installed space-heating (electrical radiator, water radiator, or floor 
heating), the distribution efficiency (ij dist ) is assumed 94%, 84%, or 
87%, respectively (Table 7). When mechanical cooling is needed, it 
will take place for a short period. Therefore, the COP cu of the nom¬ 
inal operating condition (25 °C outdoor air temperature) is used 
in presenting the annual performance. The implementation of a 
flat-plate solar thermal collector reduces the domestic hot water 
demand Qdhw by 4 Qdhw. The DHW energy saving (4 Qdhw) is sim¬ 
ulated by IDA ESBO [31] and checked by TRNSYS (Type lb) for 
different sizes of solar collector (Fig. 4). The solar thermal collector 
is orientated to south and tilted by 45°. These angles are selected to 
achieve maximum DHW energy saving. This is in accordance with 


Fig. 5. Annual electricity production [KWh] per one square meter of PV. 


[27], which indicates that the best orientation in Finland is south, 
southeast or southwest, and the optimal tilt angle is in between 30 
and 70 degrees. 

The hourly photovoltaic (PV) electricity production is simulated 
also by IDA ESBO [31 ] and checked by TRNSYS (type 194). Paramet¬ 
ric analysis (Fig. 5) shows that the maximum electricity production 
from one square meter PV module can be achieved by 50° tilt angle 
and 15° azimuth angle (from south to east). Because of the house 
orientation, the optimal azimuth is not allowed. Therefore, the PV 
collector is oriented to south. To achieve the best performance from 
south orientation, the PV is sloped by 45 degrees. The annual elec¬ 
tricity production is only 4% less than the maximum that could be 
achieved by the optimal azimuth and tilt angles. 

Fig. 6 shows the annual useful electricity generated by differ¬ 
ent photovoltaic panel sizes. The useful electricity generation is 
the production which matches hourly the total house electric¬ 
ity use, including lighting, appliances, fans, pumps, and heating. 
Fig. 6 presents three cases without flat-plate solar thermal collector 
installation. Case 1 shows the electrical demand and useful produc¬ 
tion if the reference house is heated by electrical radiators, and an 
electrical boiler is used to deliver domestic hot water (DHW). In this 
case the house electrical annual demand is 126.5 l<Wh/m 2 a. Case 2 
assumes that the house and the DHW are heated by a GSHP. Case 
3 shows that no electricity is used for heating, oil boiler or district 
heating (DH) to cover all heating demands. 


2.4.2. The difference in life-cycle cost (dLCC, second objective 
function) 

The life-cycle cost (LCC) is one of the commonly used tools to 
assess the financial viability of projects, and it has been widely 
used in energy-efficient building research [15,21-23,49,51-57]. 
LCC finds the assessment of long-term costs, as it takes into 
account discounted cash flows and covers the entire lifetime of the 
investment. From the private perspective, the LCC is the sum of 
the present value of investment and operating costs for building 
and service systems, including those related to maintenance and 
replacement, including taxes, over a specified calculation period. 
Since the current investigation aims to compare different designs in 
the specified solution space, the absolute value of the LCC is not cal¬ 
culated, but the difference (dLCQ) between the LCC for any design 
(LCC|) and that for the reference one (LCC r ). dLCC is the second 
objective function 

dLCQ = LCC/ - LCC r (2) 
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Fig. 6. Electricity demand and useful PV-electricity production where different heating systems are implemented for house with the reference envelope (Section 2.3). 
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LCC = E‘9 + Cj + MC + OC + C 


j= i 


J=5 


where IC is the investment costs of the eleven investigated design 
variables (Tables 1-8), RC is the replacement cost of the replaced 
building elements and systems (e.g., window, shading, heat recov¬ 
ery unit, etc.), MC is the maintenance costs of the heating systems 
(Table 7), IC,RC, and MC are described in Section 3.2, OC is the oper¬ 
ating cost of energy (Section 3.2), and C is a constant for other 
costs such as construction and design cost, i denotes indexes for 
the design solution, j is an index for the design parameter (e.g., 
building insulation, window type etc.). 

In this way, there is no need to include cost data for all 
components of the building but only the differences produced 
by the variation of specified parameters between the reference 
design and any other investigated design. According to the Direc¬ 
tive 2010/31/EU [8], the life-cycle cost analysis is performed for 
a 30-year calculation period. Longer calculation periods are not 
recommended, as beyond a 30-year timeframe, assumptions on 
interest rates and forecasts for energy prices become very uncer¬ 
tain [9]. It is worthwhile to mention that no disposal costs are taken 
into consideration for components with lifetimes longer than the 
calculation period. 

In line with the draft’s cost-optimal regulation [58], 3% real 
interest rate (r) and 2% energy price escalation rate (e) are used 
to demonstrate the introduced optimization method (Sections 
3.1-3.3). Sensitivity analysis with different r and e is performed to 
present the influence of the economic assumptions on the cost- 
optimal combination of building elements and heating/cooling 
system (Section 3.3). 


2.4.2A. Investment, maintenance and replacement cost calculation. 
The investment cost (IC) is the summation of the investment cost 
of wall, roof, and floor insulation, cost of window, shading cost, 
additional labour cost of improving the building tightness, the 
investment cost of the ventilation heat recovery unit, the instal¬ 
lation cost of the heating system, the cost of cooling system (if 
used), and the installation cost of solar systems. The replacement 
cost (RC) is calculated for the building elements and systems which 
have a lifespan shorter than the calculation period (30 years): win¬ 
dows, shading elements, heat recovery unit, cooling unit, solar 
thermal collector, and the inverter of the photovoltaic system. The 


replacement costs are discounted, assuming the replacement will 
take place in the middle of the calculation period. 

RQ = IC /; x(l+r)- n / 2 (3) 

where k is the subscript of the replaced building element, r is the 
real interest rate and n is the lifespan 

No maintenance costs are assumed for the replaced elements. 
Maintenance costs (MC) are assumed for the heating systems 
(Table 7), including the replacement costs of some system parts. 


2.4.22. Energy operating cost calculation. Operating costs include 
the annual fees of the electrical and the district heating (where 
used) connections and the annual heating/electrical energy costs. 
The annual connection fees (Table 6) are discounted by a discount 
factor (a), calculated from [59], considering the real interest rate (r) 
and the calculation period (n) into account 


1 - (1 + r)~ n 
r 



Different energy prices P (Table 7) are used to calculate the 
energy cost. The calculation of the electrical energy cost is based 
on the daytime/night time rates, including taxes. The daytime rate 
is applied on Mondays-Fridays between 07:00 and 20:00. The night 
time rate is applied outside these times. The energy costs are dis¬ 
counted 


OC — a e P + a e P DHWdeiivered + a e P Elejeiivered (5) 


1 - (1 + r e )~ n 


r - e 


where 

a e \ discount factor (a) takes into account the escalation rate of 
the energy price as well. 

r e : real interest rate (r) including the effect of the escalation rate 
of the energy price. 

The current study assumes there is no change in the building 
energy demand over the calculation period. 

In Finland, the energy consumption per dwelling has been rather 
steady over the years (from 1995 to 2007) when normalized for 
outdoor temperature [60]. Concerning the electric energy demand 
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of European residential users, various studies have been carried out 
[61-63]. The ODYSEE Project [63] showed that the electric demand 
per dwelling did not vary significantly in the decade between 1996 
and 2005 since, even though household appliance consumption 
declined, electric appliance penetration and use increased. 

2.5. Using a multi-stage optimization method 

In order to find the optimal trade-off between the defined 
objective functions (PEC and dLCC), multi-stage optimiza¬ 
tion method is proposed to explore more than 3xl0 9 
(16x8x13x3x3x4x3x2x4x31x71) combinations of 
the design-variable options. Considering the impact of the design 
variables on the objective functions, the exploration is performed 
in three stages: 

• Stage 1 aims to find the optimal combinations (packages) of 
the design variables which influence the thermal performance 
(heating, cooling, comfort) of the house i.e., building envelope 
parameters (Tables 1-4) and heat recovery unit (Table 5). 

• Stage 2 assesses the economical and environmental viability 
of implementing the offered primary heating/cooling systems 
(Tables 6 and 7) to the optimal combinations found in Stage 1. 

• Stage 3 investigates improving the economical and/or envi¬ 
ronmental viability of the optimal combinations of building- 
envelope parameters and HVAC systems assessed in Stage 2. 
Stage 3 addresses the RESs (Table 8 ) as supplementary systems 
for heating and/or electricity production. 

The stages provide the decision maker a trade-off relation 
between the environmental and financial viability of the studied 
building. The decision maker can choose the cost-optimal and nZEB 
solutions according to his or her preferences 


3. The optimization stages: results and discussions 


3.1. Optimal combinations of building envelope and heat 
recovery design options (Stage-1) 


The aim of the stage is to find representative energy-efficient 
building designs, disregarding the heating, cooling, and energy- 
supply systems. In order to achieve this, the space-heating energy 
demand of the house and the present worth (PW) of the influ¬ 
encing measures (insulation, building tightness, window type and 
shading option, and heat recovery type) are minimized, while 
a penalty function is applied when the summer-comfort crite¬ 
rion (DH 27 > 150°Ch) is violated. Thermal demand is minimized 
because it is the major demand in residential buildings, particu¬ 
larly in the cold-climate EU countries [64]. The PW presents the 
initial and replacement costs of the addressed measures 


PW 


= ]Tic, + ]Trc, 


i=1 


i=l 



The investment costs (IC) are shown in Tables 1 -5. The replacement 
costs (RC) are calculated by Eq. (3). According to the Finnish building 
code D3 [45], degree-hours (DH 27 ) are used to measure the summer 
overheating risk 


i=8760 

DH27 = 'y ] dT 27 At 

i=i 



dT 2 7 = (I, — 27) Tj — 27 > 0 
dT 27 = 0 when Tj - 27 < 0 


where Tj is the mean air temperature [°C] at the warmest zone 
(upper floor) and At is a one-hour time-period [h]. 

Since it is computationally expensive to explore (simulate, 
assess, then compare) all the available combinations of the 
building-envelope parameters (Tables 1-4) and the heat-recovery 
units (Table 5), automatic simulation-based optimization method 
is used. The simulation time is reduced by using a simplified three- 
zone model (Fig. 2b) for a multi-zones predefined house (Fig. 2a). 
The simplification’s details are described in a previous work [7]. To 
save time, a modified multi-objective genetic algorithm PR.GA is 
used to speed up the exploration. The algorithm is a combination of 
deterministic and a controlled elitist GA (a variant of NSGA-II [65]) 
from MATLAB toolboxes (2008b). The modified algorithm is able to 
address discrete and continuous variables, avoid repetition, keep all 
the iterations in archive and use them in non-dominated sorting 
processes. The algorithm was developed [ 66 ] and used success¬ 
fully in previous building and HVAC studies [7,67,68]. The idea of 
the algorithm is that instead of initiating the GA optimization with 
random population, the algorithm performs a preparation phase 
PR attempting to find diverse solutions close to optimums and pro¬ 
vides them to the GA as seeds. By this, the algorithm reduces the 
random behavior of the GA and converges to the Pareto-optimal 
sets faster than simple GAs would. The preparation concept can be 
adopted in different ways (e.g., parametric analysis) as shown in 
[69]. 

In the current optimization stage, the PR_GA is employed with a 
manual preparation phase (PR). The GA optimization started with 
an initial population (20 individuals), including 10 random and 
10 manually chosen diverse building-envelope and heat-recovery 
combinations: the reference design, the cheapest and the most 
expensive combinations, and seven intermediate low-energy com¬ 
binations. The intermediate low-energy combinations have lower 
space-heating demands than the cheapest solution and lower 
investments than the most expensive one. The intermediate combi¬ 
nations are chosen manually with incremental H-values in between 
the lower and upper bounds (Table 1 ). The intermediate combina¬ 
tions have the same building tightness and shading option as the 
reference design. 

Aiming to achieve high quality results, GA performed 40 
generations, considering all the design variables as discrete. 
The H-value’s boundaries of wall, roof, and floor construc¬ 
tions are divided by a constant step ( 0 . 02 W/m 2 K) to 16, 8 , 
13 options, respectively. The total building-envelope and heat- 
recovery combinations (16x8x 13x3x3x4x3 = 179712) are 
explored effectively by only 800 simulation runs. Lower num¬ 
ber of simulations was applicable. However, accurate results were 
desired. 

All the simulations are shown in Fig. 7. The optimal ones are 
classified into two groups. Group 1 satisfies the summer overheat¬ 
ing criterion (Eq. (7)), while Group 2 does not. Groups 1 and 2 
consist of 19 and 13 solutions, respectively. Group 2 solutions are 
not eliminated, as non-comfort solutions, because they could be 
feasible when the mechanical cooling option (Table 6 ) is added. 
The RES (e.g., photovoltaic) could improve the economical feasi¬ 
bility of the mechanical cooling solutions by covering a portion of 
their electricity demands. The evaluations show that investing in 
ESMs (Tables 1 -3 and 5) has a potential to reduce the space-heating 
demand from about 50 to 25 kWh/m 2 a. The influence of the ESMs 
and shading options (Table 4) on the summer overheating level 
- presented by the DH 27 - is shown in Fig. 8 . With the lowest- 
cost shading option (internal-blind shading, option 4), acceptable 
overheating levels (DH 2 7 < 150°Ch) cannot be achieved. External 
or outer panes shading (options 1 and 2 ) can achieve the comfort 
levels. However, they reduce the passive heating opportunity. From 
a thermal-comfort point-of-view, a blind between the inner panes 
(option 3) is a suitable shading option, if extreme building tightness 
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Fig. 7. Stage 1 optimization evaluation: the annual space-heating energy demands 
versus the present worth (PW) of different combinations of building-envelope 
parameters (Tables 1-4) and heat-recovery unit (Table 5). 
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Fig. 9. Financial and environmental viability of the optimal building and HVAC 
designs before and after the RES improvements (Stage-3). 


and/or insulation levels are not used. The cost-optimal building 
designs Group 1 and 2 have shading option 3 and 4, respectively. 

3.2. Combining the heating/cooling systems with optimal 
building designs (Stage-2) 

Stage 2 assesses the financial and environmental viability of 
Stage 1 cost-optimal building designs (Group 1 and 2; Fig. 7), 
addressing the offered primary heating systems (Table 7). The 
supplementary systems (e.g., solar systems) are out of the scope 
for the current stage. The stage implements mechanical cool¬ 
ing (Table 8) only for the overheated building designs (Group 2). 
Short simulation (three-month simulation) is employed to calcu¬ 
late the cooling demands during the summer season for only 13 
building solutions (members of Group 2). The heating demands 
are taken from the pre-simulated results (Stage-1). Fig. 9 shows 
the assessment results, classifying them according to the imple¬ 
mented heating/cooling systems and building group. The results 
show that although Group 2 solutions have lower space-heating 
energy demands than Group 1 ones (Fig. 7), they are not eco¬ 
nomically feasible due to their need for mechanical cooling, which 
increases not only the investment cost but also the operating one. 
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The results also show that the financial and the environmen¬ 
tal aims do not necessarily contradict each other. The highest 
capital cost heating system (GSHP) is the cost-optimal one in 
term of LCC. Annual fees, service expenses, and energy costs 
reduce, respectively, the economic feasibility of the other heating 
options: district heating (DH), oil boiler (OB), and electrical heat¬ 
ing (EH). The global cost-optimal building design requires more 
ESM implementations than the reference one. The cost-optimal 
and the reference ones have a space-heating demand of 36.8 and 
461<Wh/m 2 a, respectively. Table 9 presents the cost-optimal com¬ 
binations of building-envelope parameters and heat-recovery unit 
for each heating system. The combinations show that the cost- 
optimal implementations of space-heating ESMs (e.g., additional 
insulation) depend on the installed heating system. For economical 
operation, the highest operating-cost heating system (EH) requires 
more ESMs than the lower operating-cost ones: OB, DH, and GSHP. 

From the environmental point of view, using a highly energy- 
efficient heating system reduces the potential of the ESM 
implementations. Generally speaking, it seems that investing in 
an energy-efficient heating system has a better financial and envi¬ 
ronmental viability than investing in extremely high insulation, 
efficient window, and/or tight envelopes. An efficient heating sys¬ 
tem reduces not only the space-heating environmental burden 
but also the DHW ones. Further environmental assessments of the 
embodied energies show that the heating systems cause marginal 
impacts compared with the entire house in the pre-operating phase 
[70]. A life-cycle assessment (LCA) of GSHP systems in Europe 
showed that in most countries, particularly in countries with an 
energy-efficient electrical grid, C0 2 savings can be achieved with 
the application of shallow geothermal systems. For instance, using 
a GSHP system only for heating achieved an average of 18% savings 
in comparison to gas-furnace heating and 35% in comparison to oil- 
fired boilers [71 ]. Using GSHP as a heating option for single-family 
houses tends to increase the costs, according to the Finnish market 
statistics (Fig. 3). Our current results encourage this trend. Seeking 
for higher economical and/or environmental building-integrated- 
system designs, the next stage continues the exploration, offering 
RES options (Table 8) as supplementary systems. 

3.3. Optimal building-integrated RES-system designs: nZEB 
solutions (Stage-3) 


Space-heating energy demand [kWh/m 2 a] 

Fig. 8. Degree hours (DH 27 ) and space-heating demands of all the evaluated solu¬ 
tions. Shading options are described in Table 4. 


This stage seeks to improve the financial and/or environmental 
viability of Stage-2 building-envelope and HVAC-system designs. 
The stage investigates the viability of implementing RES options as 
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Table 9 

The reference design and the cost-optimal solutions for each heating system, presented in Fig. 9. 


Heating system 
option (Table 7) 

PEC (kWh/m 2 a) 

dLCC (€/m 2 ) 

Space heating 
(l<Wh/m 2 a) 

(/-value 

(W/m 2 K) 

Wall Roof 

Floor 

Window type 
option 
(Table 2) 

Heat recovery 
option 
(Table 3) 

Building tightness 
option (Table 4) 

Ref. design 

215 

0 

46 

0.170 

0.090 

0.170 

1 

1 

1 

EH 

178 

-33 

28 

0.095 

0.080 

0.157 

2 

3 

2 

OB 

149 

-83 

34 

0.131 

0.087 

0.171 

1 

2 

2 

DH 

112 

-57 

35 

0.154 

0.090 

0.157 

1 

2 

2 

GSHP 

103 

-104 

37 

0.142 

0.083 

0.157 

1 

2 

1 


supplementary systems for heating and/or electricity production. 
Since Stage-2 solutions (Fig. 9) have different potential for improve¬ 
ment via the RES options (Table 8), the integrated designs have 
a probability to compete with each other. The integrated designs 
are combinations of 128 (19x4 + 13x4) Stage-2 designs and 2201 
(31 x 71) RES options. The integration composes 281,728 solu¬ 
tions. In order to reduce the optimization effort, pre-investigation 
is performed to reduce the solution space by eliminating the 
unpromising solutions. The pre-investigation assessed the viability 
of adding the largest size (70 sq m) of PV panels to all Stage-2 solu¬ 
tions. The assessment showed that it is not possible to improve the 
financial and/or the environmental viabilities of the direct electric¬ 
ity (EH), oil boiler (OB), and district heating (DH) solutions further 
than the GSHP ones. The GSHP is an indirect electrical heating 
system. Therefore its solutions have a higher potential for improve¬ 
ment by the PV electricity production than the OB and DH solutions. 
Solutions with direct electricity heating system (EH) have a higher 
RES improvement potential than the ones with a GSHP. However, 
due to the efficiency of the GSHP, the EH cannot compete with it. 
Similarly, the solar thermal system shows higher economical fea¬ 
sibility in integrating with the electricity-based heating systems 
than the fuel-based ones (OB or DH). Based on the pre-investigation 
results, the OB, DH, and EH solutions are eliminated from the solu¬ 
tion space of Stage-3 exploration. 

The solution space is reduced further by eliminating the GSHP- 
dominated solutions. The pre-investigation (adding the smallest 
size 1 sqm of the thermal and photovoltaic systems, sequentially, 
to all the GSHP solutions) shows that expensive RES options cannot 
improve the economic performance of the GSHP-dominated solu¬ 
tions when compared with the cost-optimal non-dominated ones. 
Eliminating the dominated solutions, with and without a separate 
cooling option, lets the investigation focus on very energy-efficient 
building designs. Inefficient energy designs are not preferable even 
if they are economically feasible [72]. The explored solutions (non- 
dominated solutions), with and without a cooling option (fronts 2 
and 1, Fig. 9), are shown in Table 10 and Table 11, respectively. 

Tables 8 and 9 show that the optimal solutions with mechanical 
cooling could compete with the optimal solutions without cool¬ 
ing in both aspects: economical and environmental. For instance, 
although solutions without cooling (e.g., solution 8; Table 11 ) dom¬ 
inate the optimality of the solutions with cooling (e.g., solution 3; 
Table 10), the latter have a higher potential for improvement by 
RES (e.g., PV). Solution 3 (Table 10) consumes more electricity-for 
cooling-and less space-heating than solution 8 (Table 11). Since 
a portion of solution 3 electricity is for cooling during summer 
day mornings, the electricity demand has a high potential to be 
reduced by PV implementation. Consequently, solution 3 might 
compete with the optimality of solution 8 if the same PV size were 
implemented for both and there were no electricity export to the 
grid. 

In order to investigate the improvement of the cost-optimal 
building and HVAC designs (fronts 1 and 2, Fig. 9), optimization is 
performed while considering the building and HVAC design as well 


as the size of the thermal and photovoltaic solar systems (Table 8) 
as design variables. Various options are taken for each design vari¬ 
ables: (1) eighteen options (Table 10 and Table 11 solutions) for 
building and HVAC designs, (2)31 collector areas (from 0 to 30 sq m) 
for a solar thermal system, and (3) 71 module areas (from 0 to 
70 sq m) for the photovoltaic system. The total number of the three 
design variable combinations is 39, 618 (18 x 31 x 71). Instead of 
evaluating all the 39,618 combinations, the optimization performs 
two times 1200 evaluations to find the cost-optimal integrated 
solutions with and without the cooling implementation, respec¬ 
tively. Linear constraint function is used to avoid solar systems 
(thermal and PV collectors) implementation larger than the house’s 
roof area (70sqm). The optimization is performed by the mod¬ 
ified multi-criteria optimization algorithm (PR_GA) introduced in 
Section 3.1. In the preparation phase (PR), a single-objective deter¬ 
ministic algorithm (Fmincon, from MATLAB 2008b Toolbox) is used 
to minimize the first objective (PEC), considering the second objec¬ 
tive (dLCC) = 0, -50, and -100 as constraint via three sequential 
optimization runs. The evaluations of the deterministic algorithm 
are kept in the archive then sorted by the non-dominated sorting 
code in MATLAB m-file. The non-dominated solutions are used as 
an initial population for the genetic algorithm phase (GA). Stage-3 
explored the 39,618 offered solutions by only 2400 evaluations (6% 
of the offered solutions). This is a significant reduction in computa¬ 
tional effort particularly when time-consuming simulation is used 
(i.e., the simulation of the integrated design solutions consumes 
time in considering the matching between energy production and 
consumption). Further reduction was applicable. However, high 
quality results were desired. Fig. 10 presents the financial and envi¬ 
ronmental viabilities of the optimal building and HVAC designs 
(packages) before and after the RES improvements. The optimiza¬ 
tion history and the Pareto optimal sets are shown in Fig. 11a. 


o 



Fig. 10. Towards nearly zero energy house: the optimisation results (a), the opti¬ 
mal combinations of ESMs and RES (b). The ESMs implementation is presented 
by its space-heating energy saving potential as a percent of the reference design 
(46 kWh/m 2 a). 
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Table 10 

Optimal solutions with mechanical cooling (Front 2, Fig. 9), heated by GSHP. 


Index 

NPE(kWh/m 2 a) dLCC(€/m 2 ) Space heating 

(kWh/m 2 a) 

H-value (W/m 2 K) 


Window type 
option 
(Table 2) 

Heat recovery 
option 
(Table 3) 

Building tightness 
option (Table 4) 

Wall 

Roof 

Floor 

1 

99.6 

-58.9 

26 

0.069 

0.070 

0.092 

2 

3 

3 

2 

99.7 

-64.0 

26.1 

0.082 

0.073 

0.092 

2 

3 

3 

3 

99.7 

-65.0 

26.2 

0.078 

0.075 

0.103 

2 

3 

3 

4 

100.0 

-76.0 

26.7 

0.101 

0.073 

0.157 

2 

3 

3 

5 

100.5 

-77.9 

27.5 

0.114 

0.087 

0.144 

2 

3 

3 

6 

101.4 

-83.6 

28.8 

0.131 

0.080 

0.109 

2 

3 

2 

7 

102.0 

-88.5 

29.9 

0.142 

0.073 

0.171 

2 

3 

2 

8 

104.0 

-93.2 

33 

0.131 

0.083 

0.144 

1 

2 

2 

9 

108.8 

-98.2 

40.8 

0.169 

0.087 

0.171 

1 

2 

1 


Table 11 

Optimal solutions without mechanical cooling (Front 1, Fig. 9), heated by GSHP. 

Index 

NPE 

(kWh/m 2 a) 

dLCC (€/m 2 ) 

Space heating 
(kWh/m 2 a) 

DH 27 (°Ch) 

H-value (W/m 2 K) 

Window type 
option 
(Table 2) 

Heat recovery 
option 
(Table 3) 

Building tightness 
option (Table 4) 

Wall 

Roof 

Floor 

1 

95.1 

-65.4 

25.7 

148 

0.069 

0.073 

0.092 

2 

3 

3 

2 

95.2 

-72.2 

25.8 

150 

0.082 

0.073 

0.103 

2 

3 

3 

3 

95.6 

-83.6 

26.4 

148 

0.122 

0.073 

0.133 

2 

3 

3 

4 

96.8 

-88.5 

28.2 

145 

0.095 

0.080 

0.157 

2 

3 

2 

5 

97.6 

-89.2 

29.4 

147 

0.101 

0.080 

0.133 

2 

2 

2 

6 

98.0 

-93.9 

29.9 

146 

0.154 

0.083 

0.144 

2 

3 

2 

7 

98.9 

-95.6 

31.3 

143 

0.154 

0.080 

0.157 

2 

2 

2 

8 

99.7 

-98.0 

32.5 

140 

0.122 

0.080 

0.157 

1 

2 

2 

9 

102.6 

-104 

36.8 

139 

0.142 

0.083 

0.157 

1 

2 

1 
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Fig. 11. Preferable cost-optimal integrations of ESMs and RES options at different energy price escalation rates. The ESMs implementation is presented by its space-heating 
energy-saving impact as a percent of the reference design (46 kWh/m 2 a). All solutions use GSHP for heating. No cooling is implemented. 
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Fig. lib shows the design-variable combinations of the improved 
optimal solutions. The building and HVAC design is presented by 
its space-heating energy saving percentage (i.e., a building and 
HVAC design with higher energy-saving measures has lower space¬ 
heating demand than the reference design). The other design 
variables (thermal and PV solar systems) are presented by their 
areas. 

According to the Directive 2010/31/EU, the minimum LCC 
solution (global cost-optimal solution) determines the minimum 
energy performance requirements. However, a slightly higher LCC 
solution could be preferable if it reduces the PEC significantly. 
Fig. 10 shows the global and preferable cost-optimal designs 
assuming r=3% and e = 2%. The difference between the LCC of 
the cost-optimal solutions is 5€/m 2 . Based on the resulted global 
and preferable cost-optimal solutions, the minimum energy per¬ 
formance level should be 103 or 92kWh/m 2 a of primary energy, 
respectively. The cost-optimal levels are respectively 40% and 47% 
lower than the 172 l<Wh/m 2 a (372 - 1.4 x Areahouse) level, which 
has come in force by the Finnish code D3 [45]. The gap between 
the cost-optimal levels and the D3 one is larger than the EPBD 
gap threshold (15%). In term of C0 2 - eq mission the cost-optimal 
energy performance level should be in-between 27.5 and 25 kg 
C0 2 -eq/m 2 a. These values are calculated by assuming 0.459 kg 
C0 2 _ eq per 1 kWh electrical energy. 

Fig. 11 shows that the solar systems have a significant potential 
to improve the financial and environmental viability of the GSHP 
building designs. Although the exploration assumes that there is 
no export to the grid, the solar thermal and photovoltaic systems 
reduced the primary energy of the house about 40 kWh/m 2 a from 
that needed in the global cost-optimal design. Investing in a solar 
system not only reduced the primary energy but also improved 
the economic feasibility of most of the Stage-2 optimal designs 
(fronts 1 and 2). The results show that investing in solar systems can 
achieve a nearly-zero-energy building with about 70 l<Wh/m 2 a pri¬ 
mary energy consumption. However, higher investing appears not 
feasible from the economic point-of-view. Energy export, economic 
incentives, and/or cost-effective innovative measures are required 
for nearly-zero-energy buildings with a lower PEC. For instance if 
export of PV surplus energy to the grid is allowable, the nZEB level 
will be as low as 40 kWh/m 2 a. 

Concerning net-zero-energy buildings, recent publications 
[15,72] raised a question: “To what level should we decrease the 
energy use by means of energy-efficiency measures before imple¬ 
mentation of renewable energy sources?”. Our results (Fig. 11) 
show that the optimal answer is difficult to find and generalize. 
Optimization is required for answering. The economic viability of 
the ESMs and renewable systems depend on each other as well as 
on the installed heating/cooling system. The optimal implemen¬ 
tation of RES systems depends on the energy consumption of the 
building and its HVAC design. Larger photovoltaic systems are more 
economically feasible for building designs with mechanical cooling 
than for designs with proper passive-cooling measures. Mechanical 
cooling increases electricity demand; however, it allows matching 
between energy generation and consumption. In addition, it allows 
extra low space-heating solutions without the overheating risk. 
With a large size PV (>50 sq m), the financial and the environmental 
viabilities of the building design with mechanical cooling system 
become close to the viabilities of systems with proper passive¬ 
cooling measures (Fig. 10). 

Fig. 11 shows also that the financial and environmental via¬ 
bility of the photovoltaic (PV) panel is higher than the viability 
of the solar thermal system. The reason is that the PV not only 
reduces lighting and appliances’ electrical demand but also reduces 
the heating requirement. All the cost-optimal building designs are 
heated by GSHP, which uses electricity in an indirect way. The 
energy-saving potential of the PV implementation increases when 



Fig. 12. Preferable cost-optimal integrations of ESMs and RES options at different 
energy price escalation rates. The ESMs implementation is presented by its space¬ 
heating energy-saving impact as a percent of the reference design (46kWh/m 2 a). 
All solutions use GSHP for heating. No cooling is implemented. 

electrical demand increases (Fig. 5). The maximum economically 
feasible size of a solar thermal collector is 15sqm, according to 
Fig. 11. Regarding storage capacity and temperature limit consid¬ 
erations, large collector sizes (>15 sq m) cannot significantly reduce 
domestic hot water demand (Fig. 4). 

4. Sensitivity analysis 

Sensitivity analysis is performed to show the impact of energy 
price escalation rates (e) on cost-optimal results. Seeking for a 
preferable cost-optimal solution, as illustrated in Fig. 10, multi¬ 
objective optimization is required. The preferable solution is not 
the global optimal one. Therefore, a single-objective optimization 
method is not suitable for finding it. The introduced multi-stage 
multi-objective optimization method is used to find the prefer¬ 
able cost-optimal solution, at the range of e: from 2 to 15% and 
with 3% real interest rate (r). To reduce the computational effort, 
mechanical cooling is not considered as an option. The above cal¬ 
culation, with r=3% and e = 2%, shows that it is not cost-optimal 
to use mechanical cooling to avoid summer overheating risk. Pas¬ 
sive techniques (e.g., shading) is proper. This assumption reduced 
the solution space by eliminating all the building envelopes which 
suffer overheating from the exploration. Fig. 12 shows the prefer¬ 
able cost-optimal integrations of ESMs and RES options at different 
energy price escalation rates. 

The results show that GSHP is the cost-optimal heating system 
at the studied range of e. The solar thermal system cannot compete 
as a cost-optimal solution because it increases both the investment 
and replacement costs. Solar thermal collectors require replace¬ 
ment at the half of the indicated life cycle (30 years), while no 
replacement is assumed for the PV. The economic viability of invest¬ 
ing in ESMs (Tables 1-5) increases proportionally with the energy 
price escalation rate from 4% to 9%. Since the ESMs have a limited 
influence on the energy consumption, their extreme levels (e.g., 
passive D-values) do not appear as cost-optimal ones. The maxi¬ 
mum ESMs implementation has a potential for 44% space-heating 
energy saving. At e > 9%, lower ESMs are used as a cost-optimal 
solution to save 41% of the SH energy. Cost-optimal investing in 
RES (e.g., PV) depends on the space-heating demand and the esca¬ 
lation price’s escalation rate (e). Investing in 11 sqm PV-panels 
is economically feasible when a low level of space-heating ESMs 
are implemented (energy saving <20%) and low energy-price 
escalation rates (<4%) are assumed. With that level of ESMs imple¬ 
mentation, 11 sq m of PV has a potential to reduce the heating cost 
supplementing the indirect heating system (GSHP) by free electric¬ 
ity. At high rates of escalation (e > 5%), space-heating energy saving 
(>20%) is economically feasible. However, there is no potential for 
PV > 7 sq m to reduce the heating cost. Generally speaking, from the 
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economic point of view, investing in ESMs and/or the RES option is 
feasible at higher energy-price escalation rates. The optimal inte¬ 
gration of ESMs and RES options is mutually dependent and also 
depends on the economic assumption. 

It is worthwhile to mention that a small number of building and 
HVAC design options could not reach the cost-optimum and might 
lead to non-preferable environmental solution. A two-objective 
optimization algorithm is required to find the preferable cost- 
optimal solution. The global cost-optimal solution can be found 
by a single-objective algorithm. However, it does not provide the 
decision maker a preferable solution with significantly better envi¬ 
ronmental burden and slightly higher investments (5% LCC). 

5. Conclusions 

A multi-stage simulation-based optimization method is 
introduced to find cost-optimal and nearly-zero-energy building 
solutions in line with the EPBD-recast 2010. The method sig¬ 
nificantly reduces the computational effort (compared with the 
exhaustive search or trial and error methods) and supports a 
transparent cost-optimality analysis, presenting the effect of the 
different design variable combinations. The method provides a 
wide-solution-space exploration in reasonable time. 

As a case study, the method is applied to find the cost-optimal 
and nZEB energy performance levels for a single-family house in 
Finland. An enormous number of combinations (more than 3 x 10 9 ) 
of building-envelopes, heat-recovery units, cooling and heating 
options, as well as solar-thermal and photovoltaic solar systems are 
explored with a suitable number of iterations (less than 3400). The 
method speeds up the exploration, not only by reducing the num¬ 
ber of evaluations but also by using pre-simulated results instead 
of running time-consuming simulations (when possible). 

The introduced method explored the defined solution space in 
three stages. In the first stage, an efficient multi-objective genetic 
algorithm (a variant of NSGA-II [65]) is used to find cost-effective 
energy-efficient combinations of the design variants which influ¬ 
ence the thermal performance of the house: the building-envelope 
(insulation thickness of external wall, roof, and floor, window type, 
and building tightness) and the heat-recovery unit. Then the sum¬ 
mer overheating risk is checked and a mechanical cooling system 
is added to the combinations which suffer unacceptable summer 
overheating levels. The second stage assesses the primary-energy 
consumption (PEC) and the life-cycle cost (LCC) of the optimal 
combinations, addressing the offered primary heating options (e.g., 
electrical heating, oil boiler, district heating, GSHP). Then a non- 
dominated sorting code is used to determine the optimal trade-off 
between the objective functions (PEC and LCC). The third stage 
attempts to improve the financial and/or environmental viabilities 
of the optimal found building and HVAC solutions by implementing 
the optimal sizes of solar- thermal and photovoltaic systems. 

The results show that for a single-family house in cold climate 
of Finland: 

• The financial and the environmental aims do not necessarily 
contradict each other. The achieved preferable cost-optimal solu¬ 
tions, at different energy price escalation rates (from 2 to 15%) 
and a calculation period of 30 years, have primary-energy con¬ 
sumption around 47% lower than the D3-2012 standard value 
(172 kWh/m 2 a). 

• The cost-optimal implementations of ESMs and RES methods 
depend on the installed heating/cooling system and its energy- 
price escalation rate. For a cost-optimal operation, a building 
with a fully electrical heating system requires more ESMs (e.g., 
additional insulation) than a building heated by a GSHP or fuel- 
based heating systems (e.g., OB and DH). PV sizes up to 20sqm 


and 15sqm are preferable cost-optimal options for houses with 
high electrical energy demands: houses with a mechanical cool¬ 
ing option and a direct or indirect electrical heating system 
(EH and GSHP), respectively. Smaller PV sizes (up to 5sqm) are 
economically feasible for houses without a mechanical cooling 
implementation and fuel-based heating systems. Solar-thermal 
system is not a cost-optimal solution, particularly for houses 
heated by fuel-based heating systems. The solar-thermal system 
has a lower economic viability than the PV system because the lat¬ 
ter reduces the most expensive energy source (electrical energy). 
The solar-thermal system not only increases the investment cost 
but also the replacement one. The life-span of a solar-thermal col¬ 
lector is often shorter than that of a PV one. For nZEB solutions, a 
solar-thermal system with collector up to 15 sq m is economically 
feasible. 

• The cost-optimal combinations of space-heating ESMs depend 
on their influence on summer overheating. Measures which 
significantly increase summer overheating (e.g., inefficient 
shading methods) are not preferable from an economic point 
of view because they could lead to a need for mechanical 
cooling. 

• Mechanical cooling is not a cost-optimal option. It increases not 
only the investment cost but also the operating cost. Proper 
shading and building tightness as well as natural cooling via oper¬ 
able windows can eliminate the summer overheating risk. The 
economical feasibility of the mechanical cooling option becomes 
close to optimal when it is integrated with a sufficient size of PV 
for free electricity from solar panels. 

• Higher energy-price escalation rates encourage investments in 
RES and/or ESMs options. This is limited by the energy-saving 
potential of the options. For instance, increasing insulation to 
“the passive house D-values” is not a cost-optimal option because 
it has a limited influence on the space-heating energy saving. A 
lower insulation level is the cost-optimal one. 

• Currently, on-site solar systems cannot contribute as a part of a 
global cost-optimal solution because of their expensive capital 
costs. However, small PV sizes can compete with the ESMs for 
achieving energy performance levels higher than the global cost- 
optimal one, with a slight increase in the LCC. 

• From the economic point of view, it is viable to achieve nZEB 
with primary-energy consumption up to 70 kWh/m 2 a. Economi¬ 
cal and environmental incentives/credits are required to improve 
the economic feasibility of solutions towards net-zero-energy 
building. For instance if export of PV surplus energy to the grid is 
allowable, the nZEB level will be as low as 40 kWh/m 2 a. 


Exploring a huge number of design options, which include RES 
technologies, HVAC systems, and commonly used/advanced energy 
saving measures, is extremely time-consuming and arduous to 
analyse. Automatic simulation-based optimization is required to 
find the optimal combinations when having the economical and the 
environmental viabilities as objective functions and the minimum 
acceptable thermal comfort level as a constraint. A multi-stage 
optimization method can reduce the exploration effort and sup¬ 
port a fully informative and transparent analysis in line with the 
EPBD-recast 2010. 


Acknowledgments 

The first author would like to acknowledge the financial support 
of the Doctoral Programme in the Built Environment (RYM-TO). 
The second author would like to thank the Academy of Finland for 
funding his work as part of a Research Fellow position. 


202 


M. Hamdy et al. / Energy and Buildings 56 (2013) 189-203 


References 

[1] W. Eichhammer, T. Fleiter, B. Schlomann, S. Faberi, M. Fioretto, N. Piccioni, S. 
Lechtenbohmer, A. Schiiring, G. Resch, Study on the Energy Savings Potentials 
in EU Member States, Candidate Countries and EEA Countries, Final Report for 
the European Commission Directorate-General Energy and Transport, 2009. 

[2] B. Wesselink, Y. Deng, Sectoral Emission Reduction Potentials and Economic 
Costs for Climate Change (SERPEC-CC), Summary report, ECOFYS, 2009. 

[3 ] T. Boermans, J. Grozinger, Economic effects of investing in EE in buildings - the 
BEAM2 Model, Background paper for EC Workshop on Cohesion policy, 2011. 

[4] M. Thyholt, A. Hestnes, Heat supply to low-energy buildings in district heating 
areas: analyses of CO 2 emissions and electricity supply security, Energy and 
Buildings 40 (2) (2008) 131-139. 

[5] L. Gustavsson, A. Joelsson, Life cycle primary energy analysis of residential 
buildings, Energy and Buildings 42 (2) (2010) 210-220. 

[6] B. Rolfsman, CO 2 emission consequences of energy measures in buildings, 
Building and Environment 37 (12) (2002) 1421-1430. 

[7] M. Hamdy, A. Hasan, K. Siren, Applying a multi-objective optimization approach 
for design of low-emission cost-effective dwellings, Building and Environment 
46(1)(2011) 109-123. 

[8] The Directive 2010/31/EU of the European Parliament and of the Council 
of 19 May 2010 on the energy performance of buildings, Official Jour¬ 
nal of the European Union, 53, 2010, available at: http://eur-lex.europa.eu/ 
LexUriServ/LexUriServ.do?uri=OJ:L:2010:153:0013:0035:EN:PDF 

[9] T. Constantinescu, Cost optimality discussing methodology and challenges with 
the recast energy performance of building directive, The Buildings Performance 
Institute Europe (BPIE), Brussels, September 2010. 

[10] T. Boermans, K. Bettgenhauser, A. Hermelink, S. Schimschar, other Ecofys inter¬ 
national staff, Cost optimal building performance requirements, Calculation 
methodology for reporting on national energy performance requirements on 
the basis of cost optimality within the framework of the EPBD, European Coun¬ 
cil for an Energy Efficient Economy (ECEEE), 2 May 2011. 

[11] European Commission, Draft Commission Staff Working Paper, Guidelines, 
Accompanying the document, Commission Delegated Regulation of supple¬ 
menting Directive 2010/31/EU of the European Parliament and of the Council 
on the energy performance of buildings (recast) by establishing a comparative 
methodology framework for calculating cost optimal levels of minimum energy 
performance requirements for buildings and building elements, In Preparation 
of a Delegated Act in Accordance with Art 290 TFEU 6 May 2011 in Brussels, 
http://www.energy.eu/directives/2010-31 -EU.pdf 

[12] EN-15459, Energy performance of buildings, Economic evaluation procedure 
for energy systems in buildings, 2007. 

[13] ISO 15686-5, Buildings and constructed assets e service-life planning e part 5: 
life-cycle costing, 2008. 

[14] L. Georges, C. Massart, G. Moeseke, A. De Herde, Environmental and eco¬ 
nomic performance of heating systems for energy-efficient dwellings: case 
of passive and low-energy single-family houses, Energy Policy 40 (1) (2012) 
452-464. 

[15] A.J. Marszal, P. Heiselberg, Life cycle cost analysis of a multi-storey residential 
net zero energy building in Denmark, Energy 36 (2011) 56005609. 

[16] M. Elkinton, J. McGowan, J. Manwell, Wind power systems for zero net 
energy housing in the United States, Renewable Energy 34 (5) (2009) 
1270-1278. 

[17] J. Kurnitskia, A. Saari, T. Kalamees, M. Vuolle, J. Niemela, T. Tark, Cost opti¬ 
mal and nearly zero (nZEB) energy performance calculations for residential 
buildings with REHVA definition for nZEB national implementation, Energy 
and Buildings 43 (11) (2011) 3279-3288. 

[18] H. Ren, W. Gao, Economic and environmental evaluation of micro CHP systems 
with different operating modes for residential buildings in Japan, Energy and 
Buildings 42 (6) (2010) 853-861. 

[19] Rakennuslehti magazine, Helsinki, Finland, 22.01.04. 

[20] R. Holopainen, M. Hekkanen, Energy renovation saving potentials of typical 
Finnish buildings, VTT’s energy renovation technologies-project, 2006. 

[21] P. Pylsy, T. Kalema, Concepts for Low-Energy Single-Family Houses, Tampere 
University of Technology, Department of Mechanics and Design, Research 
Report 1, 2008. 

[22] A. Hasan, M. Vuolle, K. Siren, Minimisation of life cycle cost of a detached house 
using combined simulation and optimization, Building and Environment 43 
(12) (2008)2022-2034. 

[23 ] K. Alanne, A. Salo, A. Saari, S. Gustafsson, Multi-criteria evaluation of residential 
energy supply systems, Energy and Buildings 39 (12) (2007) 1218-1226. 

[24] K. Alanne, N. Soderholm, K. Siren, I. Morrison, Techno-economic assessment 
and optimization of Stirling engine micro-cogeneration systems in residential 
buildings, Energy Conversion and Management 51 (2010) 2635-2646. 

[25] A. Saari, T. Kalamees, J. Jokisalo, R. Michelsson, K. Alanne, J. Kurnitski, Finan¬ 
cial viability of energy-efficiency measures in a new detached house design in 
Finland, Applied Energy 92 (1012) (2012) 76-83. 

[26] Eurostat, Housing statistics, available at: http://epp.eurostat.ec.europa.eu/ 
statistics_explained/index.php/Housing_statistics, 2009 (accessed 02.04.12). 

[27] D5 National Building Code of Finland, Calculation of power and energy needs 
for heating of buildings, Guidelines 2012, Ministry of the Environment, Depart¬ 
ment of the Built Environment, DRAFT 28 September 2010, available in English 
at http://www.puntofocal.gov.ar/notific_otros_miembros/fin43_t.pdf 

[28] D2 National Building Code of Finland, Indoor Climate and Ventilation of 
Buildings, Regulations and Guidelines 2012, Ministry of the Environment, 


Department of Built Environment, DRAFT 28 SEPTEMBER 2010, available in 
English at http://www.puntofocal.gov.ar/notific_otros_miembros/fin41_t.pdf 

[29] J. Shemeikka, A. Laitinen, Specification of RET-single family house, Technical 
Research Centre of Finland VTT, Building and Transport, Version 2.0; 2005 (in 
Finnish). 

[30] IDA-ICE, IDA Indoor Climate and Energy, http://www.equa.se/ice/intro.html 

[31 ] IDA ESBO, http://www.equaonline.com/esbo/IDAESBOUserguide.pdf 

[32] T. Kalamees, K. Jylha, H. Tietavainen, J. Jokisalo, S. Hornets, R. Hyvonen, S. Saku, 
Development of weighting factors for climate variables for selecting the energy 
reference year according to the EN ISO 15927-4 standard, Energy and Buildings 
47 (4) (2012) 53-60. 

[33] C3 National Building Code of Finland, Thermal insulation in a building, regula¬ 
tion, Ministry of the Environment, Housing and Building Department, 2010. 

[34] Rakennustutkimus RTS Oy, available at: http://www.motiva.fi/en/building/ 
choice_of_heating_system/ (accessed 02.04.12). 

[35] RIL 249-2009, Low Energy Building Residential Buildings. 

[36] SG Isover Oy-Finland, Private communication, 2 March 2006. 

[37] Statistics Finland, Building Cost Indices, available at: http://www.stat.fi/til/rki/ 

[38] Tamar, awnings and blinds, www.tamar.fi (accessed 02.04.12). 

[39] Motiva, Specialist in Energy and Material Efficiency, http://www.motiva.fi 
(accessed 02.04.12). 

[40] NIBE Industrial Group, http://www.nibe.fi (accessed 02.04.12). 

[41] Helsingin Energia, http://helen.fi (accessed 02.04.12). 

[42] Gebwell Ltd. - Ecological Heating Solutions, http://www.gebwell.fi/eng/ 
(accessed 02.04.12). 

[43] K-RAUTA, http://www.k-rauta.fi (accessed 02.04.12). 

[44] Allergia-apu, Indoor air expert, http://www.allergia-apu.fi (accessed 02.04.12). 

[45] D3 Finnish code of building regulation, Energy efficiency of buildings, Regula¬ 
tions and guidelines 2012, Helsinki, Finland; 2011 (in Finnish). 

[46] Statistics Finland, available at: http://www.stat.fi/til/ehkh/tau_en.html 

[47] Helsingin Energia, available at: http://helen.fi/index_eng.html 

[48] European Commission (EC), Energy Technologies, Knowledge, Perception, 
Measures (Eurobarometer) 2006. 

[49] M. Leckner, R. Zmeureanu, Life cycle cost and energy analysis of a net 
zero energy house with solar combisystem, Applied Energy 88 (1) (2011) 
232-241. 

[50] Solar Buzz Solar Market research and analysis, NPD Group company, available 
at http://solarbuzz.com/facts-and-figures/retail-price-environment (accessed 
2.04.12). 

[51 ] J. Morrissey, R. Horne, Life cycle cost implications of energy efficiency measures 
in new residential buildings, Energy and Buildings 43 (4) (2011) 915-924. 

[52] A. Chel, G. Tiwari, A. Chandra, Simplified method of sizing and life cycle cost 
assessment of building integrated photovoltaic system, Energy and Buildings 
41 (11)(2009) 1172-1180. 

[53] B. Agrawal, G. Tiwar, Life cycle cost assessment of building integrated 
photovoltaic thermal (BIPVT) system, Energy and Buildings 42 (9) (2010) 
1472-1481. 

[54] G. Verbeeck, H. Hens, Life cycle optimization of extremely low energy dwellings, 
Journal of Building Physics 31 (2) (2007) 143-177. 

[55] T. Nielsen, Optimisation of buildings with respect to energy and indoor envi¬ 
ronment, PhD Thesis, Technical University of Denmark, Denmark, 2002. 

[56] M. Hamdy, A. Hasan, K. Siren, Optimum design of a house and its HVAC 
systems using simulation-based optimization, International Journal of Low- 
Carbon Technologies 5 (3) (2010) 120-124. 

[57] M. Hamdy, M. Palonen, A. Hasan, Implementation of Pareto-Archive NSGA-II 
Algorithms to a nearly Zero-Energy Building Optimization Problem, Building 
Simulation and Optimization (BS012), IBPSA-England, Loughborough Univer¬ 
sity, UK, 2012. 

[58] European Commission, Energy Directorate General, Meeting Document for 
the Expert Workshop on the Comparative Framework Methodology for Cost 
Optimal Minimum Energy Performance Requirements, In Preparation of a Dele¬ 
gated Act in Accordance with Art 290 TFEU 6 May 2011 in Brussels. 

[59] S. Fuller, S.Petersen, Life-cycle costing manual for the federal energy man¬ 
agement program, US Department of Commerce, National Institute of 
Standards and Technology, NIST handbook 135, 1995 ed., available at: 
http://fire.nist.gov/bfrlpubs/build96/PDF/b96121 .pdf 

[60] S. Elvas, L. Gynther, Energy Efficiency Policies and Measures in Finland, 
Motiva Oy, Helsinki, September 2009, available at: http://www.odyssee- 
indicators.org/publications/PDF/finland_nr.pdf 

[61] ADEME, CCE, CRES, Odense Elforsyning Net A/S, Politecnico di Milano, Servizi 
Territorio, End-use metering campaign in 400 households of the European 
Community-Assessment of the Potential Electricity Savings, Project EURECO, 
2002. 

[ 62 ] eERG, end-use Efficiency Research Group, Results of the survey campaigns elec¬ 
tricity consumption at 110 homes in Italy-Load curves of major household 
appliances and lighting, Politecnico Milano 2004. 

[63] ADEME, intelligent energy Europe. Energy Efficiency Trends and Policies in the 
EU 27, Lessons from the ODYSSEE MURE projectt, 2009. 

[64] M. Bianchia, A. De Pascalea, P. Spinab, Guidelines for residential micro-CHP 
systems design, Applied Energy 97 (2012) 673-685. 

[65] K. Deb, Multi-Objective Optimization using Evolutionary Algorithms, John 
Wiley & Sons, Chichester, UK, 2001, ISBN 0-471 -87339-X. 

[66] M. Hamdy, A. Hasan, K. Siren, Combination of optimization algorithms for a 
multi-objective building design problem, in: 11th International Building Per¬ 
formance Simulation Association IBPSA, Glasgow, UK, 2009. 


M. Hamdy et al. / Energy and Buildings 56 (2013) 189-203 


203 


[67] M. Hamdy, A. Hasan, K. Siren, Optimum design of a house and its HVAC 
systems using simulation-based optimization, International Journal of Low- 
Carbon Technologies 5 (3) (2010) 120-124. 

[68] M. Hamdy, A. Hasan, K. Siren, Impact of adaptive thermal comfort criteria 
on building energy use and cooling equipment size using a multi-objective 
optimization scheme, Energy and Buildings 43 (9) (2011) 2055-2067. 

[69] M. Turrin, P. Buelow, R. Stouffs, Design explorations of performance driven 
geometry in architectural design using parametric modeling and genetic algo¬ 
rithms, Advanced Engineering Informatics 25 (2011) 656-675. 


[70] L. Yang, R. Zmeureanu, H. Rivard, Comparison of environmental impacts of 
two residential heating systems, Building and Environment 43 (6) (2008) 
1072-1081. 

[71 ] D. Saner, R. Juraske, M. Kiibert, P. Blum, S. Hellweg, P. Bayer, Is it only CO 2 that 
matters? A life cycle perspective on shallow geothermal systems, Renewable 
and Sustainable Energy Reviews 14 (7) (2010) 1798-1813. 

[72] A.J. Marszal, P. Heiselberg, J.S. Bourrelle, E. Musall, K. Voss, I. Sartori, A. 
Napolitano, Zero energy building-a review of definitions and calculation 
methodologies, Energy and Buildings 43 (4) (2011) 971-979. 



